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ABSTRACT 

We determine simple analytical conditions for combined radiative and thermal 
equilibrium between the X-ray emitting plasma and cold reprocessor in active galactic 
nuclei (AGN) and galactic black holes (GBH). These conditions determine the pressure 
at the transition zone, which is not arbitrary in the static situation. The conditions 
for static solution are derived analytically for different forms of plasma heating and 
for cooling provided by Compton cooling, bremsstrahlung, and (optionally) advective 
cooling. 

We conclude that if Compton heating is the only heating mechanism, we always 
achieve static equilibrium between the phases. For constant volume heating and ra- 
diatively cooled plasma the static solution is never achieved and evaporation or con- 
densation takes place. However, static solutions with no evaporation/condensation are 
found for radiative cooling supplemented with lateral advection. Similar results are 
obtained for a two-temperature plasma. In the case of a general prescription of me- 
chanical heating (Q^ = hoP"^T~'^) of radiatively-cooled one-temperature plasma, we 
found that static solution is only possible, when < s < 3/2. 

Such conditions apply to all models, like disc/corona solutions, ADAF flows at 
their outer ADAF/disc boundary or cold clumps embedded in a hot medium. These 
models have to be reconsidered and supplemented either with additional condition 
for the pressure at the transition zone or with equation describing the mass exchange 
between the phases. 

Key v^rords: galaxies: active - accretion, accretion discs - black hole physics - binaries 
- X-rays; stars - galaxies:Seyfert, quasars - X-rays. 



1 INTRODUCTION 

The broad band spectra of AGN and GBH clearly show 
that the accreting matter close to a black hole consists of 
two phases. Relatively cool, optically thick phase has the 
characteristic temperature of order of 10^ K in AGN and 10^ 
K in GBH, and the electron temperature of the hot optically 
thin phase is about 10^ K (for a review, see Mushotzky, Done 
& Pounds 1993; Tanaka & Lewin 1995). 

The description of the coexistence of the hot and cold 
medium has therefore to be incorporated into all models of 
the accretion flow, independent from their assumed geome- 
try for cold and hot gas. 

Most papers on this subject model only the radiative 
equilibrium within the hot and cold medium, including the 
radiative coupling between the two phases (see e.g. Haardt & 
Maraschi 1991, Zdziarski, Lubiriski & Smith 1999). It leaves 
unnecessary arbitrariness in the modeling, e.g. the value of 



the transition radius in ADAF based models or the fraction 
of the energy liberated in the hot phase in coronal models. 

However, at the boundary between the hot and cold 
medium there is an extreme temperature gradient and the 
flow of the heat due to conduction, which is proportional to 
this gradient, is essential within the boundary zone. Even if 
the plasma coexist in radiative equilibrium, the conduction 
usually leads to a time evolution of the gas: either evapo- 
ration of the cold phase or condensation of the hot phase 
proceeds (e.g. Begelman & McKee 1990, McKee & Begel- 
man 1990). This time evolution is specifled by the efficiency 
of conduction. Efficient evaporation may lead to disappear- 
ance of the innermost part of the disc, as shown by Meyer 
& Meyer-Hofmeister (1994) for the case of accretion onto a 
white dwarf. 

The role of conduction in the context of AGN or 
GBH was studied only in a few papers so far. Maciolek- 
Niedzwiecki, Krolik & Zdziarski (1997) calculated the ver- 
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Figure 1. The heating and coohng functions niultiphed by the 
conduction flux, normahzed to their values at T = Ty^^t^ for the 
case of Compton heating. The conductive/radiative equilibrium 
means that there are three crossections between the two curves 
and the areas between the two curves above and below the inter- 
mediate crossection are equal. 



tical temperature profile for an upper part of the corona 
with direct heating of electrons proportional to the den- 
sity. Rozaiiska (1999) calculated the full temperature profile 
for a Compton heated corona. Krolik (1998) and Torricelli- 
Ciamponi & Courvoisier (1999) studied the case of a cloud 
model. Krolik (1998) considered the conductive balance of 
the clouds with constant heating per unit volume. Torricelli- 
Ciamponi & Courvoisier (1999) calculated numerically the 
temperature profile across the cold cloud/hot inter-cloud 
medium for a few special cases of the heating and for 
coohng function including bremsstrahlung as well as other 
atomic processes important at lower temperature. Finally, 
Mineshige et al. (1998) discussed the quiescent state of X-ray 
novae taking into account the disc evaporation and DuUe- 
mond (1999) considered the transition zone in the case of 
one-temperature, advection-cooled hot corona above an ac- 
cretion disc. 

In this paper we rediscuss the problem of condition im- 
posed by the effect of conduction onto the boundary between 
the cold and hot phase. We show in a simple analytical way 
which forms of heating lead to well defined condition for 
pressure at the transition zone, supplementing the condi- 
tion of purely radial equilibrium deeply within the hot layers. 
Heating mechanisms which lead to heating constant per unit 
volume (like the aP scaling of Shakura & Sunyaev 1973), or 
give too steep drop of the heating with the temperature, do 
not allow for any static boundary between the two phases 
and in such models explicit time evolution (e.g. evaporation 
rate of the cold phase) should be included. 



2 GENERAL CONDITIONS OF 

RADIATIVE/CONDUCTIVE EQUILIBRIUM 

The hot plasma and cold gas may coexist both in hydro- 
static and radiative equilibrium if they are cooled by two 
different mechanisms supporting different temperature equi- 
librium. However, this equilibrium in general is perturbed by 
the presence of conduction and either the cold phase would 
evaporate or the hot phase would condensate due to the heat 
exchange by electrons. Additional condition has to be satis- 
fied to sustain the two-phase equilibrium, and generally this 
condition can be formulated as a condition for the pressure 
at the hot/cold plasma transition zone. 

The conductive heat flux Fcond is proportional to the 
temperature gradient. Therefore, it is extremely impor- 
tant in all those boundary regions, where the temperature 
changes by three or more orders of magnitude over the small 
distance determined by the Field length (Field 1965). We 
consider here in details the case of plane-parallel geometry 
of the transition zone but appropriate analysis can be re- 
peated in other geometries (see Appendix A) 

The classical expression for Fcond is given by: 

i^ rri5/2dT 

tcond = -KqJ —^, (1) 

dz 
where the temperature T changes in the direction of z coor- 
dinate. For a plasma of cosmic abundance, the conductivity 
constant is equal to (Draine & Giuliani 1984) : 



K.0 — 5.6 X 10 4>c 



(2) 



where (j)c is the factor corresponding to reduction in the 
mean free path due to magnetic fields and turbulence {(j)c = 
1 for equal ion and electron temperature). 

This classical expression has to be replaced with the 
expression for saturated conduction fiux if the mean free 
path of electrons becomes larger than the characteristic scale 
high for the temperature. However, in the context of AGN 
or GBH this is not necessary (e.g. Rozanska 1999). 

The energy balance across the boundary between the 
hot and cold medium in the static case is given by: 

^^ cond 



dz 



Q^ 



A, 



(3) 



where Q^ determines the (radiative or non-radiative) heat- 
ing of the gas and A describes the radiative cooling. The 
boundary conditions are set by the requirement that the 
conductive fiux vanishes both deeply in the cold zone and 
deeply in the hot zone: 



^ cond\Zcold) U ^cond\Zhot)^ 

and the hot zone remains in radiative equilibrium: 

Q'^{zhot) = A{zhot)- 



(4) 



(5) 



Equation (H) with boundary conditions (M) and n se- 
cure that neither evaporation nor condensation take place 
and the boundary between the two phases remains at its 
equilibrium position. 

Condition (kl) can be equivalently expressed in the in- 
tegral form: 

((5+ - A)dz = 0. (6) 

old 

Those conditions can be implemented numerically into 
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the complete set of equations determining the transition 
zone between the hot and cold medium. However, we can 
obtain simplified analytical conditions determining the equi- 
librium transition zone. 

In most cases studied in detail in the context of strat- 
ification in AGN or GBH the rapid change of the temper- 
ature (and density) by a few orders of magnitude happens 
practically under the condition of almost constant pressure, 
as noticed by Torricelli-Ciamponi & Courvoisier (1998) and 
Krolik (1998), and used in their approach. 

The condition of constant pressure allows to relate the 
density of the gas to its temperature; 

_ PuniH 



Tk 



(7) 



with the value of molecular weight ^ = 0.5 for cosmic chem- 
ical composition. Boltzmann constant is denoted as k and 
the mass of hydrogen atom as niH- 

This gives the opportunity to express all quantities as 
functions of the temperature if the heating term also de- 
pends entirely on the density, temperature and constant fac- 
tors not varying across the boundary region. Equation (pi) 
can be solved in general form by the change of variables from 
ztoT (and from d/dz to FcondT'^^^ /tiod/dT): 

l-Thot 

FlUTcou) - Fl„^{Thot) = 2 / (Q+ - A)Acor"/'dr.(8) 

This formal solution can be directly used to formulate 
the conditions of radiative/conductive equilibrium. We fur- 
ther simplify the conditions noticing after Krolik (1998) that 
the cold medium temperature is much lower than the hot 
medium temperature, so it can be approximated by zero. 
Imposing boundary conditions (M) and (Bl) we then finally 
arrive at two conditions: 



(Q+ - A)K.oT^''^dT = 



and: 



Q^{Thot)^A{Thot). 



(9) 



(10) 



Those two conditions determine both the hot medium 
temperature and the pressure at the transition zone if the 
heating prescription is assumed and the cooling mechanisms 
specified. 



3 SPECIAL SOLUTIONS 

The cooling mechanisms acting in AGN and GBH are rel- 
atively well known. Here we consider Compton cooling and 
bremsstrahlung: 



Ak 



A{p,T) = Fradi^.s ^pT + Bp^T'', 



(11) 



and expressing the density p by the pressure P constant 
through the transition zone we have: 



A[p,l ) = t radices -^ —F+B( -—) F 1 , 



rrieC^ 



(12) 



where k^s is opacity coefficient for electron scattering , rrie - 
mass of electron, c - velocity of light, and B - bremsstrahlung 
cooling constant, which equals B = 6.6 x lO^^erg/scm^g^. 
Since we approximated the cold temperature by zero we 



do not have to include complex bound-free and line cool- 
ing important at the lowest temperatures. As was noticed 
by Torricelli-Ciamponi & Courvoisier (1998), the detailed 
description of the processes at lower temperatures is not es- 
sential since the conductive flux is proportional to T^i'^ and 
the integrand is strongly weighted towards That- 

We also neglect here the synchrotron cooling. However, 
if the magnetic field in the hot plasma is important the syn- 
chrotron cooling term should be included. We discuss the 



special case of advection cooling in Section 3.3 



As for heating, the physical process leading to observed 
high temperature of electrons is not known. Several param- 
eterization of this process are used in modeling. Therefore 
we study here the most frequently used assumptions consid- 
ering both one-temperature and two-temperature plasma. 
In this last case we have to supplement the equations with 
electron/ion Coulomb coupling. 

3.1 Compton heating 

Compton heated plasma typically achieves temperature of 
order of 10^ — 10* K, depending on the spectral shape of the 
incident radiation. This is not enough for X-ray source in 
AGNs and GBHs (see introduction), but Compton heated 
plasma is often considered, e.g. as an outer corona in AGN 
or as confining medium for broad line clouds in AGN, so it 
is interesting to check how equations H and lid work in this 
simplest example. 

The term describing Compton heating can be written 
in the usual form: 



^ — ^ radl^e 



4k 



■pTic, 



where Tic is the Inverse Compton temperature. 
Equations (BI) and ( |lO| ) lead to the conditions: 

7 



Tho 



and 



15 



Tic, 



15 "^^ ''" BprnHtneC^ "^' 



tT 



3/2 



(13) 



(14) 



(15) 



These conditions (see also Figure hh reproduce well the 
requirements for the coexistence of the hot and cold phase 
based purely on radiative stability, although the coefficients 
are slightly difi^erent. We can introduce the ionization pa- 
rameter H of Krolik, McKee & Tarter (1981): 

I — I ^ rad 



cF 



(16) 



The requirement for the pressure leads to the same scaling 
with the Inverse Compton temperature 



H' <xT, 



-3/2 



(17) 



as determined by Begelman, McKee & Shields (1983). 



3.2 Constant volume heating and a P heating 

Looking for additional heating mechanism, which can in- 
crease the temperature of an X-ray source up to the lO^K 
we consider aF heating of the hot plasma (Shakura & Sun- 
yaev (1973), Witt, Czerny & Zycki (1997), Krolik (1998)). 
Since we assume that the transition zone is narrow enough 
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to be approximated by P = const both heating parameteri- 
zations are equivalent. 

Denoting the constant heating as Qtonat^ conditions (bl) 
and (hol) are satisfied only if: 

(18) 



and 

That = OO, Phot = 0. 



(19) 



Finite required value of the pressure is accompanied by 
condition that bremsstrahlung cooling is negligible. It effec- 
tively means that any finite temperature will lead to the vio- 
lation of one of the two conditions. Radiative equilibrium far 
from the cold medium would lead to unbalanced conductive 
flux and subsequent evaporation of the cold medium. Static 
solutions are never achieved. 

If only the integral condition (i.e. bl) is used we repro- 
duce the results expressed by Krolik (1998) as his H„et = 
7/4 (his equation 6) or by Witt, Czerny & Zycki (1997) as 



S oc T^ 



-3/2 



Such a condition for pressure minimizes the 



effect of evaporation, therefore its use is not entirely mean- 
ingless, but it does not provide strict static equilibrium. 

Similar results were obtained numerically for the case of 
spherically symmetric cloud by Torricelli-Ciamponi & Cour- 
voisier (1998). They used better description of the cooling, 
extending it down to the temperatures when other atomic 
processes are more efficient than bremsstrahlung. Their so- 
lutions derived for a number of heating functions satisfied 
the integral condition but displayed strong second derivative 
of the temperature, d?T /dz^ oc Q^ — A (see their Appendix 
A). 

3.3 Constant volume heating Avith 
radiative/advective cooling 

The models with radial advection and vertical stratification 
of the hot/cold accreting gas allow to parameterize the ad- 
vection cooling term as 






(20) 



with dimensionless coefficient S replacing the radial deriva- 
tives of thermodynamical parameters of the flow and (here 
constant) Tvir being the local virial temperature (see e.g. 
Abramowicz et al. 1995, Janiuk, Zycki & Czerny 1999). 

Both equations (tel) and (hol) can now be satisfied. The 
condition for the temperature can be written in the implicit 
form 



Thot _ 27 Q+-Ac 
T 

-* VI'. 



35 g^ 



(21) 



where Ac is the Compton cooling term (i.e. first term of 
Equation Pl2|) proportional to the pressure P. The value of 
the pressure is given by 



8k\Q+-Ac)TZ 



3/2 



(22) 



35/i2m|^B 

If the Compton cooling is negligible, those conditions can 
be used directly. Otherwise, explicit solutions for That and 
P have to be derived numerically. Evaporation of the disc 
underlying the advective flow is not expected. 



DuUemond in his paper (1990) considered hot corona 
on the top of the cold disc with exactly the same heating 
and cooling processes like in this section, but for arbitrary 
value of the pressure expressed by his dimensionless constant 
C. He concluded that disc has to evaporate, because the 
heat flux deep inside cold phase stays non zero (Eq. 32 of 
DuUemond 1990). His calculations were done for C — 0.15, 
but we can find that for C = 0.1548 both conditions are 
satisfied: energy balance equation and zero heat flux inside 
the disc, so the evaporation is not required. 



3.4 General mechanical heating and radiative 
cooling 

Mechanical heating of the hot phase can in most cases be 
parameterized in the general form as: 



= hoP'^T- 



(23) 



(see Torricelli-Ciamponi & Courvoisier 1998 and the refer- 
ences therein). 

Conditions (bl) and (|lO|) reduce to: 



rpS 

J- hot 



7hom^c^ P" 



3/2 



12^mH Fradfies 7/2 - S ' 



and 
P 



„_2 _ 3-B/i rriu 7/2 - S yS-| 

hot 



(24) 



(25) 



4/iofc2 s 

Combining those two equations we find prescription for pres- 
sure: 



plijt 



SBfi^m^ f ThoT 



4:hok'^ yl2 yLmnFradl^e 

3/2 -s / 7/2- 
s V3/2- 



(i-A) 



(26) 



All equations above can be satisfied, leading to positive 
values for both the pressure P and the hot plasma temper- 
ature, Thot only if 

< s < I (27) 

We see that the problem met by Torricelli-Ciamponi 
& Courvoisier (1998) to satisfy both the integral condition 
across the cold/hot boundary and the radiative equilibrium 
in the hot phase arose because of unfortunate choice of the 
index a (equal 1.5 and 2.25 in their case, see also Appendix 
A). 



3.5 Two-temperature plasma 

Efficient predominant heating of ions may lead to a two- 
temperature plasma, with electron temperature Tg lower 
than ion temperature, Ti. The heat exchange of plasma is 
provided by the Coulomb coupling and possibly other mech- 
anisms like plasma waves. For simplicity, we can assume ex- 
clusive heating of ions. Coulomb coupling as the only heat 
transfer mechanism, and we can neglect the contribution to 
the pressure from the electrons. 

For any ion heating parameterization Q^{Ti) we have 
following heating and cooling functions for electrons: 
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0+ = |^V(^,-T.)T-^/^ 

2 niH 

4fcT 

A(p, Te) = Fradl^esP ^ + +5(:; 



Ti , r>.2rnl/2 



-^ i;i7 



+ Bp-'T^ 



(28) 



(29) 



where D is Coulomb coupling constant equal D — 2.44 x 
lO'^^lnA and InA ~ 20. The density can be determined by 
the pressure: 



PiriH 



(30) 



which allows to express all quantities as functions of Te. The 
equilibrium conditions are given by Equations (bl) and (|lO|), 
as before, only for the electron temperature. These equations 
have to be solved numerically, after specifying Q'^ {Ti) 

We considered in detail a few special cases. 

For Q+ = const, without advection, the result is the 
same as in the case of one-temperature plasma (see Section 
3.2). If radiation balance is satisfied the integral condition is 
never satisfied for a finite hot medium temperature. Strictly 
static solutions do not exists. 

In the case of Q"*" oc p, without advection, the solution 
was found for a representative set of parameters, again simi- 
larly to the case of one-temperature plasma (see Section 3.1 
Solutions could also be found for Q"*" — const, but with the 
advection term included. 



4 DISCUSSION 

Static solution for two-phase equilibrium in AGN and GBH 
with radiative/conductive energy exchange between phases 
may, or may not exist, depending on assumed heating and 
cooling mechanisms operating in the hot phase and in the 
transition zone. 

In the case of the Compton heated corona (Rozahska 
1999) we always obtain static solution, and the required 
pressure depends on Inverse Compton temperature. 

Compton heating, however, is not sufficient to heat up 
the X-ray emitting plasma in those objects and other mech- 
anisms must operate, extracting the gravitational energy of 
the accreting material. The exact path of the energy flow is 
not known so we considered typical heating parameteriza- 
tions. 

We found that the evaporation/condensation equilib- 
rium (i.e. static solution) for radiatively cooled plasma is 
possible if the heating is a moderately decreasing function 
of the temperature (see Eqs. M and M). 

Frequently used assumption that the heating is pro- 
portional to aP (Shakura & Sunyaev 1973) does not lead 
to equilibrium for a radiatively cooled plasma. Accretion 
disc covered by such a corona would evaporate continuously. 
However, if coronal advection is allowed, equilibrium may be 
restored. 

Detailed computations of such an equilibrium, partic- 
ularly for a two-temperature plasma, have to be performed 
numerically (see Section 3.5), so the existence of solutions for 



any particular set of parameters should be checked. A num- 
ber of very interesting models of the accretion fiow have 
been recently proposed, e.g. Esin, McClintock & Narayan 
(1997), CoUin-Souffrin et al. (1996), Poutanen, Krolik & 
Ryde (1997), Krolik (1998), Nayakshin, Rappaport & Melia 
(1999), and we can apply the approach outlined in the 



present paper to these models. Our analysis is not applica- 
ble when coronal heating occurs in small compact regions, 
like in case of magnetic field reconnection. 

In realistic two-phase fiow we probably have two mech- 
anisms leading to the mass exchange between the hot and 
the cold phases. First mechanism should be directly related 
to the dynamics of the flow, like heating the plasma in mag- 
netic loops or instabilities within the cold disc, leading to 
disruption of the cold flow and hot optically thin plasma 
production. This mechanism would be responsible for the 
appearance of the hot plasma in the flrst place. The second 
mechanism operates when the hot plasma already exists: 
evaporation/condensation due to conduction leads to fur- 
ther mass exchange. This second mechanism was the sub- 
ject of the present paper and we formulated the conditions 
under which no further mass exchange proceeds. 

However, in real situation the study of the appropriate 
time scales based on adopted dynamical model should be 
done. As an example, we consider a scenario with energy 
in both phases released via viscosity dissipation (situation 
calculated in Section 3.2). We can estimate accreting and 
isobaric cooling time scale defined as follows: 



Vr 



^cool 



5 p 



(31) 



(32) 



when r is the distance from the central object, Vr - radial 
velocity. To compare those two time scales we assume n^A 
to be Compton cooling and we use Eq. hq to determine con- 
stant volume heating. The ratio of both time scales can be 
expressed by Keplerian velocity and the velocity of sound 
as: 






'-' -^ vir 

5^r2" 



(33) 



''Cool 

For accreting corona, if the temperature is much smaller 
than the virial temperature, this ratio is always much big- 
ger than 1 and radiative/conductive equilibrium should be 
achieved. Since in the case of constant heating without ad- 
vection we did not find a static solution it means that rapid 
evaporation or condensation will result in transition to be 
an one phase medium. However, if the corona temperature 
is close to the virial temperature, the evaporation time scale 
is comparable to the viscous time scale, so the two phase 
medium may exist but continuous evaporation of the cold 
material should be calculated and included into the model. 
Similar discussion of the thermal and viscous time scale 
has to be repeated for other two-phase accretion models in 
order to check whether the initially non-equilibrium condi- 
tions will be rapidly changed into practically static solu- 
tions by adjusting the pressure and the temperature pro- 
file through evaporation/condensation and radiative heat- 
ing/cooling. 
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APPENDIX A: APPLICATION FOR THE 
SPHERICALLY SYMMETRIC SYSTEM 

Similar analysis can be repeated for a spherically symmetric, 
dense cold cloud embedded in the hot, optically thin plasma, 
as it was done by McKee & Begelman (1990) and Torricelli- 
Ciamponi & Courvoisier (1999). In such case the existence 
of static solution strongly depends not only on the form of 
the heating function, but also on the size of a cloud, Rd. 

For spherically symmetric case the energy balance 
across the boundary between hot and cold medium in the 
static case is given by: 



dr 



A, 



(Al) 



where r is the radial coordinate equaling zero at the 
cloud center. Adopting classical expression for Fcond = 
—kqT^ dT/dr and assuming the same boundary conditions 
as in plane-parallel case, the equilibrium between the cloud 
and the surrounding medium is described by: 

° r\Q+ - A)/s:oT'/^dr = 0. (A2) 

Now, Tcoid is the temperature of the cold cloud. Static solu- 
tion also requires that far from the transition interface the 
radiative equilibrium is achieved. 

Those two conditions: Eq.A2 and Eq. |lO| cannot be 
solved in general since the first condition requires the knowl- 
edge of T{r) explicitly in the spherical geometry. Never- 
theless, we can analyze the situation qualitatively, study- 



ing asymptotic behavior of clouds with radii very small or 
very large in comparison to the Field length Xp, defined as 
the characteristic maximum length over which the energy 
exchange by thermal conduction is negligible and radiative 
processes are dominated (Field 1956). 

It was shown by McKee & Begelman (1990) for the case 
of radiatively heated clouds, and by Torricelli-Ciamponi & 
Courvoisier (1999) for more general heating functions, that 
when Rci « \f the static solution for spherically symmet- 
ric cloud is never achieved, such a small cloud should always 
evaporate. But McKee & Begelman (1990) found that for 
Rci ^ Xf there is the value of pressure at which the mass 
exchange is stopped and static solution condition can be 
satisfied. For higher pressure a cloud will condense (conden- 
sation occurs as a negative sign of evaporation rate). The 
range of pressures for which a cloud can condense is wider 
for cloud radii much larger than the Field length. 

Those results in the limit of small cloud radii 
Rci/Xf « 1 were confirmed analytically by Torricelli- 
Ciamponi & Courvoisier (1999) who made appropriate Tay- 
lor expansion and showed that static solution is never 
achieved. 

We made a similar Taylor expansion in the opposite 
limit, for \f /Rd « 1- For such purpose the expansion of 
r is done around point Tl which is the second crossection 
of cooling and heating functions (see Fig. Ill) . For the value 
of temperature Tl the radius is denoted as Rd, so the ex- 
pansion is given by: 



r{T) = R^i + JL,(^T~T1) 



Rti+^R'd^{T-Tl) 



Substituting the expansion to the Eq. 

''°*(Q+-A)«;oT''"rfT+ 



(A3) 



A2 wc obtain: 



+ ^^ / (0+-A)(T-Tl)^oT^/^rfr = 0. (A4) 



4 dr 
RTidT 

The second term is always positive so the conditions 
derived for a plane parallel geometry are replaced by in- 
equalities in the case of a spherical geometry. Quantita- 
tive estimate of this correction can be done using the con- 
cept of the Field length. We show an example of the re- 
sult for the case of mechanical heating m — 1 and s — 1 
(see Sec. 3). Since the Field length is inversely proportional 
to the characteristic temperature scalelength, the derivative 
dlnr/d lnT — Xf / Rd- We estimate that the second integral 
in Eq. A4 is of the order of Q+kqT^^^- Now this condition 



together with fact that for T^ot the radiative equilibrium 
is achieved, allow to determine the temperature of the hot 
phase and the pressure of the transition. Comparing to the 
temperature and pressure in case of the plane parallel ge- 
ometry when Rd = oo, we obtain: 

.\f. 



ThotiRd) = ThotiRd = oo)(l - 10 



Rc: 



(A5) 



and 

P(Rd)=P{Rd=^)il + ^^^). (A6) 

In the limit of \f /Rd << 1 for spherical clouds the static 
solution is achieved, if the solution exists in the case of a 
plane parallel geometry. 
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Between those two limits the situation is more compli- 
cated. To find the value of pressure together with spherical 
cloud radius for which the static solution is satisfied the 
proper numerical computations should be done and we will 
address this problem in our future work. 

This paper has been processed by the authors using the 
Blackwell Scientific Publications M^njX style file. 
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